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PREFACE 


This  report  is  submitted  in  compliance  with  the  requirements  of 
Contract  F30602 - 76 -C -00 54  and  covers  work  carried  out  during  the  period 
15  October  1975  to  15  September  1976.  Related  work  under  a previous 
contract  ( F30602  - 75  - C -00  1 2 ) is  reported  in  RADC-TR-75  - 185  (July  1976) 
and  RADC-TR-76- 189  (June  1976). 
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of  RADC  where  the  meteorological  and  profile  data  was  reduced.  Copies 
of  the  radiosonde  data  were  provided  by  RADC.  We  thank  J.  Ochs  and 
R.  Lawrence  of  the  NOAA  Environmental  Research  Laboratories  for  pro- 
viding us  with  the  results  of  their  airborne  experiments  and  numerical 
modeling.  H.  Yura  of  Aerospace,  Inc.  , was  very  cooperative  in  providing 
us  with  the  results  of  his  theoretical  calculations.  A.  MacGovern  and 
J.  LeFebvre  of  Itek,  Inc.  , were  particularly  helpful  with  regard  to  the  RTAM. 
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INTRODUCTION 


1.  BACKGROUND  AND  OBJECTIVES 

Atmospheric  turbulence  has  a significant  impact  on  the  operation  of 
large  aperture  optical  systems.  Degradations  observed  include  loss  of 
resolution  when  imaging,  increased  difficulty  in  detecting  dim  objects,  loss 
of  tracking  accuracy,  as  well  as  a number  of  other  effects.  The  evaluation 
of  the  performance  of  any  system  and  the  specification  of  operating  para- 
meters for  systems  designed  to  eliminate  or  minimize  these  degradations 
require  detailed  empirical  information  about  the  turbulent  environment  in 
which  they  must  operate. 

This  program,  Turbulence  Environment  Characterization,  has  as  its 
objective  the  collection,  processing,  analysis  and  interpretation  of  data  from 
a variety  of  experiments  in  order  to  provide  a characterization  of  the  tur- 
bulent environment.  The  site  of  the  measurements  is  the  ARPA  Maui  Optical 
Station  (AMOS)  atop  Haleakala  on  the  island  of  Maui,  Hawaii.  Experimental 
instrumentation  which  has  been  or  will  be  employed  include  the  Seeing 
Monitor,  the  Real  Time  Atmospheric  Measuring  System  (RTAM),  the  Star 
Sensor,  a differential  Hartmann  sensor,  a small  aperture  photometer, 
microthermal  probes  and  an  acoustic  sounder.  During  a previous  program, 
Turbulence  Characterization  and  Control,  much  of  this  instrumentation  was 
installed  and  several  periods  of  data  collection  were  carried  out.  Details 
of  this  work  can  be  found  in  Refs.  1 and  2. 

During  the  period  covered  by  this  report,  a considerable  amount  of 
data  was  collected  and  analyzed.  In  addition,  several  upgrades  to  the  in- 
strumentation were  accomplished  and  an  initial  series  of  tests  with  the 
RTAM  were  completed.  The  overall  program  status,  accomplishments  and 
the  most  significant  results  and  conclusions  are  given  below.  More  detailed 
discussions  are  given  in  the  balance  of  the  report.  Section  II  summarizes 
the  experimental  operations  and  the  status  of  all  instrumentation.  Sections 
III,  IV  and  V summarize  the  Seeing  Monitor,  Star  Sensor  and  meteorologi- 
cal data,  respectively.  The  complete  and  detailed  data  is  given  in  the 


1.  M.  G.  Miller  and  P.  F.  Kellen,  Turbulence  Characterization  and  Control. 
Interim  Technical  Report,  Contract  #F30602- 75- C- 001  2 (Avco  Everett 
Research  Laboratory),  Rome  Air  Development  Center  Technical  Report 
#RADC-TR-75-  185  (July  1975),  (A015759)  . 

2.  M.  G.  Miller,  P.  L.  Zieske  and  G.  Dryden,  Turbulence  Characterization 
and  Control,  Final  Technical  Report,  Contract  # F30t>02-  75-  C-  00  1 2 (Avco 
Everett  Research  Laboratory’),  Rome  Air  Development  Center  Technical 
Report  #RADC-TR-  76-  189  (June  1976),  (A027155). 
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Appendices.  Section  VI  discusses  various  data  comparisons  and  includes 
comparisons  with  theory  and  other  experimental  results.  The  final  section 
gives  a brief  analysis  of  some  of  the  data  collected  during  the  initial  period 
of  RTAM  testing. 

2.  PROGRAM  STATUS  AND  ACCOMPLISHMENTS 

As  of  the  date  of  this  report,  the  status  of  the  instrumentation  is  as 
follows . 

• Data  collection  with  the  meteorological  instrumentation,  Seeing 
Monitor  and  Star  Sensor  has  been  reduced  to  a level  such  that  the 
normal  observatory  personnel  can  operate  these  systems  as 
scheduled. 

• The  PDP-8  data  processing  system  has  been  operating  reliably. 
Input/output  capability  will  be  upgraded. 

• Because  of  operational  failures,  the  acoustic  sounder  was  returned 
to  RADC  for  complete  evaluation.  It  is  expected  back  early  in 
November  1976. 

• The  new  tape  recorder  has  been  installed  and  evaluated. 

• The  RTAM  has  undergone  an  initial  period  of  AMOS  testing. 

• Operations  with  the  differential  Hartmann  sensor  and  small  aper- 
ture photometer  have  been  discontinued. 

Important  accomplishments  include  the  following: 

• An  intensive  period  of  data  collection  during  late  November  and 
early  December  1975. 

• Establishment  of  routine  data  collection  activities  starting  in 
April  1976. 

• A series  of  comparative  measurements  between  the  Star  Sensor 
and  small  aperture  photometer. 

• A series  of  noise  evaluation  tests  with  the  Star  Sensor. 

3.  RESULTS  AND  CONCLUSIONS 

Two  hundred  and  twenty- eight  ten  minute  samples  of  rQ  have  been 
obtained  on  twenty- four  nights  with  the  Seeing  Monitor  yielding  the  following 
information. 

• The  ensemble  has  an  approximate  Gaussian  distribution  with  a mean 
of  9.6  cm,  a standard  deviation  of  2.2  cm  and  a range  of  (5.3  - 17.8) 
cm  at  5000  A. 


• The  statistical  significance  of  a given  value  is  estimated  to  be 
better  than  ± 2%  with  a precision  of  better  than  ± l/2%.  Estimated 
relative  uncertainty  due  to  known  errors  is  (5  - 10)%. 

"•  A potentially  major  uncertainty  is  the  absolute  calibration  of  the 
entire  system.  No  simultaneous  and  independent  measure  of  rQ 
has  yet  been  obtained. 

• A variety  of  behavior  over  periods  of  hours  has  been  seen,  including 
near  constant  conditions,  large  variability  over  time  periods  of 
minutes  and  trends  towards  both  improving  and  degrading  seeing. 

• On  some  occasions,  substantial  atmospheric  non- stationarity  is 
indicated . 

Star  Sensor  data  was  collected  on  30  occasions.  One  hundred  and 
twenty- six  upper  atmospheric  profiles  and  218  twenty- minute  samples  of 
the  log- amplitude  variance  were  collected.  These  data  yield  the  following 
information. 

• Significant  changes  in  the  profiles  are  seen  in  time  scales  of  20  to 
40  minutes.  Profiles  tend  to  flatten  and  in  some  cases  indicate 
increased  levels  of  turbulence  in  the  vicinity  of  the  meteorological 
tropopause . 

• Atmospheric  non- stationarity  effects  were  sufficiently  strong  to 
invalidate  approximately  40%  of  the  profiles. 

• The  aperture  average  (36  cm)  twenty  minute  log- amplitude  variances 
have  a mean  of  5.8  x 10'^  with  a range  of  (1.5-  28)  x 10~4.  The 
ensemble  of  twenty- four  second  variances  collected  on  a single  night 
show  an  approximate  Gaussian  distribution. 

• Simultaneous  Star  Sensor  and  small  aperture  photometer  data  show 
agreement  in  log-  amplitude  variances  of  approximately  0.  5 x 10'^ 
for  the  most  directly  comparable  results. 

• The  results  of  the  noise  evaluation  tests  have  lead  to  a more  com- 
plete characterization  of  noise  properties  then  previously  existed. 
The  operational  limit  for  reliable  information  has  been  established 
as  third  to  fourth  stellar  magnitude.  The  effect  of  noise  on  the 
reduced  profile  has  not  yet  been  determined. 

Data  comparisons  lead  to  the  following  results. 

• Agreement  between  the  empirical  profile  and  the  Seeting  Monitor 
data  is  not  good.  Spatial  variability  of  the  turbulent  field  could  be 
the  cause  of  this  discrepancy. 

• Some  correlation  is  seen  in  the  Seeing  Monitor  and  Star  Sensor  data, 
particularly  relative  to  atmospheric  non- stationarity . However, 
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because  these  two  instruments  tend  to  emphasize  different  portions 
of  the  atmosphere,  a strong  correlation  should  not  be  anticipated. 

Theoretical  modeling  of  the  profile  based  on  gross  meteorological 
data,  yields  rough  agreement  with  the  data. 

The  agreement  of  this  data  for  both  scintillation  and  correlation 
scale  agrees  quite  reasonably  with  other  reported  data. 

The  initial  tests  of  the  RTAM  yield  the  following  conclusions: 

Current  performance  levels  would  provide  data  which  would  be 
extremely  difficult  to  process  for  high  frequency  MTF  information. 

It  should  be  possible  to  implement  RTAM  for  an  independent  mea- 
sure of  r Q for  comparison  with  the  Seeing  Monitor. 

Recause  of  several  effects,  it  is  unlikely  that  useful  high-frequency 
phase  information  can  be  obtained. 
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II.  EXPERIMENTAL  OPERATIONS  AND  INSTRUMENTATION 


1.  OPERATIONS 

During  the  contract  period,  several  types  of  operations  were  carried 
out.  The  first  of  these  was  an  intensive  data  collection  period  which  started 
in  early  November  and  continued  until  early  December  1975.  The  measure- 
ment systems  used  during  this  period  were  the  microthermal  probes,  acous- 
tic sounder,  Star  Sensor,  Seeing  Monitor  and  the  b-29  photometer  equipped 
with  a small  aperture.  The  microthermal  probes,  Seeing  Monitor  and 
photometer  outputs  were  processed  using  the  PDP- 8 data  processing  sys- 
tem. The  analog  signals  from  these  sensors  were  sampled  1 350  times 
during  approximately  8.25  minutes  followed  by  a data  printout  which  re- 
quired approximately  1.75  minutes.  This  resulted  in  a basic  10  minute 
data  processing  cycle.  More  detailed  discussions  of  these  various  systems 
can  be  found  in  Ref.  2.  The  Seeing  Monitor  and  photometer  were  mounted 
on  the  b-  37  and  b-  29  telescopes,  respectively.  Because  these  two  tele- 
scopes share  a common  mount,  data  from  these  systems  could  not  be  col- 
lected simultaneously.  Hence,  alternate  10  minute  cycles  were  used  to 
obtain  data  with  these  two  devices.  The  twenty  minute  Star  Sensor  cycle 
was  synchronized  with  the  other  systems  yielding  a single  data  cycle  for 
every  two  PDP-8  cycles.  Typical  data  runs  lasted  for  three  hours  and  were 
started  at  various  times  during  the  night.  Successful  data  runs  with  some 
or  all  of  these  systems  were  accomplished  on  ten  occasions.  During  the 
week  of  17-21  November,  personnel  from  the  NOAA  Environmental  Research 
Laboratories  collected  data  with  a light  aircraft  instrumented  with  a micro- 
thermal  probe.  Coupled  with  other  instrumentation,  this  yields  an  empir- 
ical estimate  of  the  turbulence  profile  above  the  site  which  can  then  be 
compared  with  the  direct  optical  measurements.  The  results  obtained 
during  this  period  are  discussed  in  other  sections  of  this  report. 

Subsequent  to  the  above  activities,  a routine  data  collection  pro- 
cedure was  established  using  the  Star  Sensor  and  Seeing  Monitor.  Data 
runs  were  typically  two  hours  in  length  and  were  implemented  using  the 
same  basic  10  minute  collection  cycle  discussed  above.  Operations  were 
generally  scheduled  on  two  nights  a week  but  were  often  preempted  by 
higher  priority  operations  or  bad  weather.  Data  was  collected  on  an  addi- 
tional twenty-one  nights  during  the  period  16  April  to  13  July  1976.  Near 
the  end  of  this  period,  atmospheric  data  using  the  technique  of  speckle 
interferometry  was  also  collected  by  a different  experimental  group.  Some 
simultaneous  data  using  the  Seeing  Monitor  was  obtained.  When  reduced, 
this  data  (speckle)  may  provide  an  independent  measure  of  r0  for  com- 
parison with  the  Seeing  Monitor  results. 

A third  type  of  activity  carried  out  during  this  period  were  special 
operations  whose  objective  was  to  establish  proper  operation  of  instrumen- 
tation, collection  of  simultaneous  and  comparative  data  and  noise  evaluation 
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tests.  Systems  involved  included  the  Star  Sensor,  b- 29  photometer  and 
the  acoustic  sounder.  Details  of  this  work  are  discussed  in  Sections  II  2 
IV  . 4,  and  IV  . 5. 

2.  INSTRUMENTATION 


The  Seeing  Monitor,  Star  Sensor,  meteorological  sensor  and  the 
PDP-8  data  processing  system  have  been  established  as  routine  systems 
and  their  operation  can  be  carried  out  by  the  normal  observatory  personnel 
on  an  as  scheduled  basis.  Maintenance,  while  time  consuming  on  occasion 
has  not  been  a major  source  of  difficulty.  Several  modifications  and  addi-  ’ 
tions  to  some  of  the  systems  have  been  carried  out  or  are  being  planned. 
Minor  software  changes  have  been  made  to  the  Star  Sensor  to  facilitate 
output  presentation  and  remote  fine  drive  controls  have  been  installed  so 
that  optimization  of  this  instrument's  operation  can  be  carried  out  without 
entering  the  dome.  A problem  with  the  Seeing  Monitor  which  has  become 
more  severe  is  the  unbalance  in  the  two  seeing  angle  outputs.  The  source 
of  this  appears  to  be  difficulties  with  the  optical  alignment  of  the  dove  prism. 
Because  studies  of  potential  non- isotropic  effects  is  not  a high  priority  item, 
the  prism  will  be  removed  before  mounting  on  the  1 .6  meter  telescope.  It 
is  expected  that  this  will  eliminate  much  of  the  problem  and  thus  allow  other 
effects  such  as  wavelength  scaling  to  be  studied.  The  rountine  meteorologi- 
cal sensor  outputs  have  been  connected  (in  parallel  to  the  PDP-8  A/D  con- 
verter) to  a digital  meter  equipped  with  a sequencing  switch  and  scaling 
electronics  located  in  the  observatory  control  room.  This  allows  a direct 
visual  reading  of  the  gross  meteorological  conditions  (wind  direction  and 
speed,  temperature  and  dew  point)  without  resorting  to  data  reduction  with 
the  PDP-8.  It  is  also  expected  that  an  automatic,  remotely  read  barametric 
sensor  will  be  added  in  the  near  future.  The  PDP-8  has  been  operating 
reliably.  However,  the  existing  teletype  input/output  has  proved  restrictive 
for  certain  types  of  work.  Therefore,  a high  speed  paper  tape  reader/punch 
will  be  obtained  and  installed. 


The  major  difficulty  encountered  has  been  with  the  acoustic  sounder. 
Although  it  appeared  that  proper  operation  had  been  established  prior  to  the 
November- December  experiments  and  raw  data  was  collected  and  recorded 
during  the  entire  period,  the  data  processing  yielded  results  which  were 
dominated  by  noise.  Subsequently,  a variety  of  tests,  modifications,  data 
collection  attempts  and  other  remedial  actions  were  taken,  none  of  which 
yielded  usable  data.  Hence  the  instrument  was  returned  to  RADC  for  a 
complete  checkout  and  evaluation.  It  is  now  apparently  operating  correctly 
and  will  be  returned  to  AMOS  in  early  November  1976. 

The  b- 29  photometer  equipped  with  a one- inch  effective  aperture, 
while  operating  correctly,  provided  very  marginal  data  during  the  November- 
December  experiments.  The  reason  for  this  was  the  excessive  amount  of 
noise  seen  in  the  data  (approximately  90%  of  the  total  measured  variance) 
without  the  capability  of  a sufficiently  accurate  measurement  of  this  effect 
The  causes  were  probably  the  high  electrical  bandwidth  ( * 4 kHz)  used  and 
a low  efficiency  of  the  entire  system.  Consequently,  only  a very  rough 


estimate  of  the  atmospheric  log- amplitude  variance  was  obtained.  Later 
in  the  year,  the  photometer  was  reconfigured  to  a 36  cm  effective  aperture 
and  scintillation  data  was  collected  simultaneously  with  the  Star  Sensor. 

The  results  of  these  experiments  are  given  in  Section  IV  . 4 . Operations  with 
this -system  have  been  discontinued. 

The  new  Bell  and  Howell  VR3700B  tape  recorder  has  been  received 
and  is  now  operational.  This  machine  has  a fourteen  channel  record  (one 
direct  and  13  FM)  and  five  channel  reproduce  (one  2 MHz  direct,  two  80  kHz 
FM  and  two  400  kHz  FM'  capability  with  one  inch  magnetic  tape.  Nine  speeds 
are  available  covering  the  range  from  15/16  to  120  ips.  This  instrument 
provides  a data  recording  and  reproduce  compatibility  between  AMOS,  the 
Everett  Laboratory  and  other  facilities  such  as  RADC.  Initial  applications 
will  be  for  recording  acoustic  sounder  raw  data  for  on-site  processing  at  a 
later  time  and  full  bandwidth  (1  kHz)  recording  of  the  Seeting  Monitor  angular 
outputs  for  processing  at  Everett. 

The  initial  AMOS  testing  of  the  Real  Time  Atmospheric  Measure- 
ment System  (RTAM)  have  been  completed.  A summary  of  these  tests  and 
an  analysis  of  some  of  the  data  is  given  in  Section  VII. 

The  portable  differential  Hartmann  device  has  not  been  used  to 
collect  data  since  August  1974.  Because  of  higher  priority  tasks  and  the 
availability  of  other  instrumentation,  no  activity  with  this  system  is  antici- 
pated at  this  time. 
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III.  SEEING  MONITOR  DATA 


1.  GENERAL 

The  Seeing  Monitor  is  basically  a va riable- spatial- frequency 
spinning-  reticle  photometer  mounted  in  the  focal  plane  of  a large  telescope 
imaging  a distant  star  through  turbulence  (Figure  1).  Assuming  a reticle 
intensity  transmission  (T)  which  is  a spatially- varying  sine  wave,  the 
output  of  a photomultiplier  (PM)  placed  behind  the  reticle  will  include  a 
term  whose  envelope  is  proportional  to  the  Fourier  transform  of  the  star 
image  (|l(f)|)  at  the  spatial  frequency  (f)  of  the  reticle.  Properly  nor- 
malized, this  quantity  is  just  the  Modulation  Transfer  Function  (MTF)  of 
the  combined  atmosphe re- telescope  optical  system  which  is  related  to  the 
Mutual  Coherence  Function.  By  utilizing  a variable  spatial- frequency 
reticle,  the  entire  MTF  can  be  scanned  as  the  reticle  rotates.  In  actual 
practice,  a logarithmic-varying,  square-wave  pattern  is  used,  however, 
the  basic  effect  remains  the  same.  A more  complete  description  of  this’ 
instrument  is  given  in  Refs.  2 and  3. 

While  a number  of  outputs  are  available  from  this  device,  the  one 
used  in  this  experiment  is  a single  parameter  characterization  of  the  width 
of  the  MTF  or  equivalently,  the  size  of  the  image.  The  development  of  this 
signal  is  indicated  in  the  top  portion  of  Figure  2.  The  PM  output  is  an  os- 
cillating signal  whose  modulation  decreases  with  time  (increasing  spatial 
frequency).  At  the  beginning  of  each  scan,  the  modulation  is  maximum 
because  the  reticle  bars  are  larger  than  the  image.  As  the  scan  continues, 
modulation  decreases  as  a result  of  the  image  being  larger  than  the  reticle 
bars.  The  signal  is  demodulated  resulting  in  an  envelope  which  is  pro- 
portional to  the  MTF.  This  signal  decreases  as  a function  of  scan  time. 
Analog  electronics  are  used  to  determine  the  time  (i.e.,  spatial- frequency ) 
at  which  the  MTF  drops  to  one-half  its  initial  (zero  spatial  frequency)  value. 
A voltage  proportional  to  this  half-MTF  frequency  is  developed  and  provided 
as  an  output  from  the  instrument.  This  voltage  changes  with  each  succes- 
sive scan  of  the  reticle  indicating  the  changing  size  of  the  image  as  a result 
of  the  variable  atmosphere.  The  system  is  duplicated  for  an  orthogonal 
channel.  Data  processing  is  shown  in  the  lower  portion  of  the  figure.  While 
tape-recording  of  the  raw  data  will  be  utilized  in  the  future  for  detailed 
temporal  characterization,  at  present  the  outputs  are  fed  directly  to  a real 


3.  C.  R.  Giuliano,  et  al,  Space  Object  Imaging.  Final  Technical  Report, 
Contract  # F30602- 74- C- 0227  (Hughes  Research  Laboratory),  Rome 

Air  Development  Center  Technical  Report  *RADC- TR-  76-  54  (March  1976). 
(A023497). 

D.  Fried,  J.  Opt.  Soc . Am.  56,  1 372  (1966). 
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Fipuro  1.  Seeing  Monitor  - Basic  Concept.  Components  shown  are  the 
Stellar  Source  (S),  Telescope  (T),  Reticle  (R)  and  Photo 
Multiplier  (P). 
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time  processing  system.  This  system  includes  an  analog  to  digital  con- 
verter, PDP-8/I  computer  and  teletype  output.  The  PDP- 8 i6  programmed 
to  calculate  the  mean  and  variance  of  the  input  signals.  Total  sample  size 
and  sampling  rate  are  variable  and  control  the  total  cycle  time  of  the  sys- 
tem. The  average  voltage  calculated  by  the  processing  system  is  then  fit 
to  a' theoretical  model  to  provide  an  estimate  of  the  correlation  scale  (rQ). 

The  theoretical  model  used  to  obtain  an  estimate  of  rQ  is  indicated  in 
Figure  3.  The  data  processing  system  averages  a voltage  which  is  propor- 
tional to  the  spatial- frequency  at  which  the  MTF  decreases  to  one-half  its 
initial  value.  The  model  used  is  the  Short  Exposure  Average  MTF.  (4)  This 
model  is  appropriate  because  the  Seeing  Monitor  automatically  eliminates 
overall  image  wander.  The  parameter  o is  assumed  to  be  one  (near  field 
approximation)  because  of  the  large  aperture  of  the  telescope.  However, 
the  model  deals  with  the  averaged  MTF  whereas  the  processing  actually 
averages  a signal  which  is  proportional  to  the  logarithm  of  the  spatial  fre- 
quency at  which  the  MTF  on  a realization  by  realization  basis  drops  to  half 
its  initial  value.  The  assumption  is  made  that  the  spatial- frequency  (<fi/2>) 
corresponding  to  the  average  voltage  approximately  equals  the  spatial  fre- 
quency (fD)  at  which  the  short  exposure  average  MTF  decreases  to  half  its 
zero  frequency  value. 

Table  1 summarizes  the  most  significant  characteristics  of  the  ex- 
perimental and  data  processing  system.  These  conditions  yield  an  estimate 
of  rD  for  each  10  minutes  of  operating  time.  While  the  data  is  collected  at 
a central  wavelength  of  6200  A,  the  results  have  been  scaled  with  wave- 
length to  5000  R. . 

2.  ERROR  ESTIMATES 

A variety  of  effects  are  present  in  the  Seeing  Monitor  which  contri- 
bute to  errors  and/or  uncertainties  in  the  reduced  data.  These  include  finite 
sample  size,  quantization,  quantum  fluctuations,  internal  noise  sources, 
telescope  MTF  and  voltage  averaging.  Another  important  aspect  of  the  data 
reduction  is  the  absolute  calibration  of  the  angular  size  of  the  image.  Ideallv, 
error  estimates  should  be  associated  with  each  of  these  sources. 

In  the  processing,  a random  voltage  is  averaged  which  potentially 
has  contributions  not  only  from  the  atmosphere  but  also  the  noise  sources. 

A basic  question, is  the  level  of  statistical  significance  associated  with  the 
finite  sample  size.  Assuming  independent  samples,  theoryl5'  yields  a 
variance  of  the  sample  mean  { - see  Figure  2), 


5.  A.  Papoulis,  Probability,  Random  Variables  and  Stochastic  Processes 
(McGraw-Hill,  N.  Y.  1965),  245. 

6.  D.  Middleton,  An  Introduction  to  Statistical  Communication  Theory, 
(McGraw-Hill,  N.  Y.  I960),  426. 
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Figure  3.  Seeing  Monitor  - Theoretical  Model.  The  upper  figure  approxi- 
mates the  Seeing  Monitor  operation  and  the  lower  one  indicates 
the  Theoretical  Model  assuming  near-field  conditions. 
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SYSTEM  CHARACTERISTICS 
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where  o is  the  true  statistical  variance  and  N is  the  sample  size  (1  350). 

The  sample  variance,  v,  is  an  estimate  of  o^.  For  the  data  reported  here, 
v was  almost  always  less  than  1 (volts)^  and  the  sample  mean  was  typically 
in  the  range  (4-6)  volts.  Hence,  the  statistical  significance  of  the  average 
voltage  (one  sigma)  is  better  than±l%.  Tracing  this  level  of  fluctuation 
through  the  instrumental  calibration  and  theoretical  model*4)  yields  an  es- 
timated statistical  significance  in  rQ  of  better  than  ± 2%.  In  many  cases  it 
was  better  than  ±1%. 

Instrumental  noise  sources  and  quantum  fluctuations  become  signifi- 
cant because  of  two  effects.  One  is  the  degradation  in  statistical  significance 
due  to  the  addition  of  a large  noise  variance  to  the  total  variance  which  is  not 
adequately  reduced  by  the  sample  size.  The  data  supports  the  conclusion 
that  this  does  not  occur.  In  addition,  a theoretical  calculation^)  for  quantum 
fluctuation  indicates  a single  realization  (one  millisecond)  root  variance  to 
mean  ratio  of  better  than  1%  for  a second  magnitude  star.  This  value  is 
well  below  the  actual  ratio  seen  in  the  data.  The  other  effect  is  an  addition 
to  the  average  voltage  (i.  e.  , bias)  which  is  not  due  to  the  atmosphere.  As 
far  as  is  known,  this  effect  does  not  occur. 

A noise  source  which  is  not  averaged  by  the  sample  size  is  the  quan- 
tization associated  with  the  data  processing.  This  is  really  not  noise,  but 
rather  the  precision  of  the  measurement.  The  ten  bit  resolution  of  the 
PDP-8  yields  a precision  in  the  voltage  output  of  approximately  0.01  volts. 
This  in  turn,  for  typical  conditions,  yields  a precision  in  the  r estimate 
of  better  than  ± l/2%. 

More  important  than  the  above  effects  are  unknowns  which  change 
the  calibration  of  the  system.  One  of  the  more  important  of  these  is  the 
telescope  MTF  ( tq  in  Figure  3).  In  fitting  the  data  to  theory,  an  aberration- 
free,  clear  aperture  of  1 . 22  M was  assumed.  This  is  obviously  not  the  case. 
Because  of  the  low  frequencies  involved,  the  central  obscuration  of  the  tele- 
scope does  not  contribute.  Aberrations  may  be  important.  Hartmann 
testing  of  the  telescope  indicates  an  80%  encircled  energy  diameter  of 
(0.7-  1.3)  arc  sec  with  a scatter  in  the  data  of  approximately  0.4  arc  sec  . 

This  information  by  itself  is  not  sufficient  to  evaluate  the  effect  on  the 
estimate  of  rQ  because  the  detailed  shape  of  the  MTF  is  not  known.  Over 
the  '■ui.ire  experimental  period,  the  maximum  averaged  one- half  MTF 
spatial- frequency  observed  was  approximately  1 (arc  sec)'*,  This  in  turn 
implies  that  the  telescope  MTF  had  to  be  at  least  0.5  at  this  spatial- frequency . 
Fitting  this  value  to  a Gaussian  function  and  calculating  the  80%  encircled 
energy  diameter  yields  a value  of  0.6  arc  sec.  While  this  value  is  low,  it 
is  not  inconsistent  with  the  Hartmann  data.  Using  the  same  Gaussian  func- 
tion to  estimate  rQ  leads  to  the  result  that  aberrations  are  not  important 
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(<  10%)  for  values  lower  than  13  cm.  However,  for  values  of  rQ  in  excess 
of  1 5 cm,  aberrations  could  result  in  estimates  substantially  lower  than 
actual.  Hence  these  large  values  should  probably  only  be  viewed  as  lower 
bounds . 


The  theoretical  model  used  to  fit  the  data  corresponds  to  the  spatial- 
frequency  at  which  the  averaged  MTF  with  overall  tilt  removed  is  reduced 
to  one- half  of  its  initial  value.  In  practice,  the  voltage  averaged  is  log- 
arithmly  related  to  the  spatial- frequency  at  which  single,  short  exposure 
MTF's  decrease  to  0.5.  Clearly,  these  two  (theory  and  practice)  are 
different.  Depending  on  the  type  of  statistics  involved  and  parameter  values, 
the  difference  could  be  substantial.  If  the  voltage  statistics  are  assumed 
Gaussian  then 
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where  (v)  is  the  average  voltage,  ov^  is  the  voltage  variance,  <C  f ] /2  is 
the  average  frequency  and  a is  the  frequency  variance.  The  numerical 
constants  arise  from  the  instrumental  calibration.  For  most  of  the  data 
reported,  the  voltage  variance  was  small  enough  so  that  its  effect  on  the 
average  frequency  was  less  than  10%;  typically,  in  the  range  of  (1-3)%.  A 
similar  effect  occurs  in  the  estimate  of  rQ. 

Evaluating  the  effect  of  frequency  averaging  is,  in  general,  com- 
plicated. It  is  related  to  the  zero  crossing  problem  which  is  common  in 
statistical  theory.  (&)  However,  because  of  the  low- frequencies  involved,  a 
simple  first-order  analysis  can  be  developed.  At  these  low  frequencies,  the 
MTF  can  be  characterized  by  the  sum  of  the  averaged  short  exposure  MTF 
and  a small,  mean-zero  fluctuation.  (?)  Because  of  the  small  variances  in- 
volved, the  averaged  short- exposure  MTF  at  the  frequency  detected  by  the 
device  (fj/2  *n  Figure  3)  can  be  expanded  in  Taylor  series  about  the  fre- 
quency at  which  it  drops  to  one- half  its  initial  value  (f  in  Figure  3).  This 
yields 


7.  Actually,  the  random  fluctuation  is  more  appropriately  modeled  as 
having  a small,  but  non-zero  mean  (see  Ref.  8).  However,  this  should 
not  significantly  effect  the  first  order  analysis. 

8.  D.  Korff,  G.  Dryden  and  M.  G.  Miller,  Optics  Comm.  _5,  1 87  (1972). 


22 


T(fl/2)!  = 0 • 5 - <rSE(f1/2)>+  (fl/2-f0)-|r<rSE(f1/2)> 


+ 2(fl/2'fo)  df2  <'tSE  (f0^  + 


'a'  / 

+ A(fl/2> 


terms  through  second^ rde r ^ d 1' ° ° ' Averaglng  this  rcsult  and  maintaining 


fo  ' <fl/2>  + 


<(fl/2-fo>2>  df2<TSE(f)> 

2 i-; — ; — r + 

df  ^tSE  ,fo» 


rrying  out  the  indicated  averages  and  derivatives  yields 


fo*  <fl/2> 


1/2>{1  - H/41  0f2/O|/2>2  j 


Assuming  the  Gaussian  model  of  Ea  (3)  vieMc:  , • 

with  typical  values  of  less  than  2%  max,mum  effect  „f  10% 

cal,br.,t™mTherse  itVude'de'fccus"  aU°  "T  u"”t.ai„,ies  in  the  system 

these  are  considered  to  be  significant.  tUrbulence-  However,  none  of 

degraded  ^addtUon^"? ‘b*  *‘”‘-dph".ri« n?*' 

internal  inconsistency  UnfoVtunate^  indicated  an 

While  some  potentially  comparative  speckle  Lll  e*iSt' 

has  not  yet  been  reduced.  Once  the  rSaM  becomes  Spe^ionaf.  TmV 


also  be  possible  to  obtain  true  simultaneous  and  independent  estimates  of 
the  half-MTF  frequency  under  identical  conditions.  However,  thi'^  will  not 
provide  a complete  system  (including  telescope)  calibration.  Data  collection 
against  planets  could  provide  such  a calibration,  at  least  for  the  large  image 
6ize~s  that  would  be  involved. 

In  summary,  the  statistical  significance  of  a given  value  of  rc  is 
estimated  to  be  better  than  ± 2%  with  a precision  of  better  than  ± 1/2%. 

Known  and  estimated  error  sources  yield  a maximum  uncertainty  in  rQ  of 
10%  with  5%  being  more  typical.  Because  of  telescope  aberrations,  large 
values  of  rQ  (>15  cm)  may  only  represent  a lower  bound.  A potentially 
major  uncertainty  in  the  results  is  the  absolute  calibration  of  the  entire 
system.  However,  existing  and  anticipated  data  may  be  able  to  provide  an 
independent  measure  of  the  correlation  scale. 

3.  RESULTS  AND  CONCLUSIONS 

During  the  period  11  November  1975  to  1 3 July  1976,  Seeing  Monitor 
data  was  collected  on  twenty-four  nights.  These  data  are  summarized  in 
Appendix  A.  Table  A- 1 gives  the  dates,  collection  period,  number  of  points 
and  the  mean  and  range  of  rQ  at  5000  R for  each  data  run.  The  various 
graphs  in  the  Appendix  give  the  temporal  behavior  of  the  data  runs.  In  total, 
228  data  points  were  collected. 

A variety  of  temporal  behavior  are  seen  in  the  results.  This  ranges 
from  almost  a constant  value  of  rQ  over  a three  hour  period  to  large  fluc- 
tuations between  points  separated  by  twenty  minutes.  Extremes  of  this  type 
of  behavior  are  shown  in  Figure  4.  The  data  of  6 December  1975  is  essen- 
tially constant  (i.e.,  less  than  5%  variation  about  the  mean)  while  that  of 
8 December  1975  shows  large  fluctuation  (almost  ±30%).  More  typically, 
variations  over  a given  data  run  are  less  than  ± 20%. 

Another  type  of  behavior  seen  were  trends  towards  improving  (in- 
creasing rQ)  and  degrading  (decreasing  rD)  conditions.  Examples  of  this 
are  shown  in  Figure  5.  On  8 July  1976,  except  for  several  points,  seeing 
conditions  degraded  in  essentially  a monotonic  fashion  from  a high  of 
ll.o  cm  near  the  beginning  of  the  run  to  a low  of  8.  1 cm  near  the  end  of  the 
run,  a variation  of  (+  21%,  - 1 6 %)  about  the  mean  value.  Just  the  opposite 
effect  occurred  on  12  July  1976.  During  this  three  hour  run,  rG  changed 
from  a low  of  7.6  cm  near  the  beginning  of  the  run  to  14.4  cm  near  the  end 
of  the  run,  a variation  of  (-  32%,  + 30%)  about  the  mean  value.  Again,  the 
trend  (towards  improved  seeing)  is  reasonably  monotonic. 

The  total  data  sample  is  summarized  in  Figure  6.  The  mean  value 
of  the  228  member  ensemble  is  9. 6 cm  with  a standard  deviation  of  2.  2 cm . 
The  range  of  measurements  were  from  5.  3 to  17.8  cm.  Thirteen  percent 
of  the  values  lie  below  7.5  cm  and  eighteen  percent  are  above  11.5  cm . 

The  straight  line  in  the  figure  is  a theoretical  Gaussian  distribution  with  the 
empirical  mean  and  variance.  As  can  be  seen,  the  fit  i6  quite  reasonable. 


CUMULATIVE  PROBABILITY  ( % ) 


Figure  6.  Seeing  Monitor  Data  - Correlation  Scale  Distribution.  The 

straight  line  is  for  a Gaussian  distribution  with  the  empirical 
mean  and  variance. 
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Another  interesting  aspect  of  the  data  is  the  level  of  fluctuation  seen 
in  the  output  voltage  during  a given  ten  minute  collection  cycle.  If  the  single 
frame  MTF  is  modeled  as  the  sum  of  a low-frequency  deterministic  portion 
(i.  e.  ..averaged  short  exposure  MTF)  and  a Rayleigh  distributed  random 
fluctuation  which  is  band  limited  by  the  telescope  aperture,  (8)  then  the 
variance  of  the  MTF  is  controlled  by  the  latter  term  and  is  proportional  to 
(rQ  D)^.  Tracing  this  through  the  data  processing  yields  the  conclusion!*?)  that 
if  the  atmosphere  is  truly  stationary  during  a collection  cycle  (i.  e.  , constant 
r0),  the  output  voltage  should  have  an  atmospheric  induced  variance  of  no 
more  than  0.  1 (volts)2.  In  the  data,  variances  of  this  size  occurred  less 
than  10%  of  the  time.  Approximately  40%  of  the  samples  have  variances 
between  (0.  1 - 0.  2)  and  another  30%  have  variances  between  (0.  2-0.  3). 

The  average  variance  for  the  entire  data  set  is  approximately  0.  25  (volts)2. 

For  the  data  of  Figures  5 and  6,  the  run  averaged  variances  are  0.  2 for 
6 December,  1.4  for  8 December  and  0.  15  for  8 and  12  July.  The  large  value 
of  8 December  implies  that  large  fluctuations  were  also  occurring  within  a 
ten  minute  collection  period.  While  some  portion  of  the  measured  variance 
is  probably  due  to  noise,  theory  and  the  observed  minimum  value  imply  that 
the  noise  contribution  is  no  more  than  0.  05  (volts)^.  Hence  the  data  supports 
the  conclusion  that  atmospheric  non- stationarity  effects  occur  much  of  the 
time  with  a characteristic  scale  of  less  than  ten  minutes.  Similar  conclusions 
have  been  reported  elsewhere.  (10) 


9.  Gaussian  voltage  statistics  are  assumed. 

10  D.  Karo  and  A.  Schneiderman,  J.  Opt.  Soc.  Am.  £>6,  1065A  (1970) 
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IV.  STAR  SENSOR  DATA 


1.  GENERAL 

The  Star  Sensor  is  a variable  spatial-frequency  photometer  which 
scans  the  scintillation  pattern  in  the  aperture  of  a 36cm  telescope  and  pro- 
vides a measure  of  the  aperture  averaged  log-amplitude  variance  and  the 
log-amplitude  variance  associated  with  a specific  spatial-frequency  which 
is  varied  as  a function  of  time.  The  assumption  of  a theoretical  model 
allows  these  data  to  be  interpreted  in  terms  of  the  strength  of  the  turbulent 
fluctuations,  Cn^,  at  seven  different  levels  in  the  upper  atmosphere.  In 
addition,  the  aperture -ave raged  log -amplitude  variance  is  available  as  an 
output.  A complete  data  collection  cycle  requires  approximately  twenty 
minutes  yielding  a reduced  profile  at  the  end  of  this  time  period.  A more 
complete  description  of  this  device  is  given  in  Ref.  11  and  12. 

During  the  period  11  November  1975  to  13  July  1976,  operations  with 
this  instrument  were  carried  out  on  thirty-one  nights.  Valid  profiles  were 
obtained  on  twenty-seven  occasions  and  scintillation  data  only  on  an  additional 
three.  A summary  of  these  data  are  given  in  Appendix  B.  Included  are  the 
average  and  range  of  scintillation,  the  complete  set  of  valid  turbulent  profiles 
and  graphs  of  the  nightly  averaged  profiles. 

2.  TURBULENT  PROFILES 

As  shown  in  Table  B-l,  a total  of  229  turbulent  profiles  have  been 
obtained.  Of  these,  24  are  not  valid  because  of  mechanical  malfunctions  or 
saturation  of  the  photomultipliers,  leaving  a total  of  205  for  which  the  instru- 
ment was  apparently  operating  correctly.  The  data  processing  is  such  that 
each  24  sec  an  average  value  of  the  aperture  averaged  log -amplitude  variance 
is  obtained  and  recorded.  During  the  twenty-minute  cycle  required  to  carry 
out  the  complete  spatial  frequency  scan,  approximately  forty  such  values 
are  accumulated.  A final  step  in  the  processing  (before  determination  of 
the  profile)  is  to  calculate  the  mean  and  standard  deviation  to  mean  ratio 
for  this  forty  member  ensemble.  This  latter  quantity  is  a measure  of 
several  effects.  Two  of  the  more  important  are  non-stationarity  of  the 
turbulent-statistics  and  the  existence  of  cloud  obscuration  during  a protion 
of  the  data  run.  Both  of  these  effects  have  serious  consequences  in  the 


11.  G.  R.  Ochs,  R.S.  Lawrence,  T.  Wang  and  P.  Zieske,  SPIE  Proceedings 
7_5,  48  (1976). 

12.  G.  R.  Ochs,  T.  Wang,  R.S.  Lawrence  and  S.  F.  Clifford,  Applied  Optics 
2^,  2504  (1976). 
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profile  determination  because  they  can  cause  the  data  collection  conditions 
to  change  substantially  during  a cycle  and  violate  one  of  the  basic  assump- 
tions of  the  theoretical  model.  Hence  a large  value  of  the  standard  deviation 
to  mean  ratio  indicates  a profile  of  questionable  significance. 

For  operational  purposes,  a value  of  this  ratio  in  excess  of  0.  2 has 
been  established  as  unacceptably  hich.  Based  on  this  criteria,  79  profiles 
were  eliminated  leaving  a total  of  126  collected  during  the  period  on  twenty- 
seven  nights.  These  are  tabulated  in  Table  B-2  and  include  the  data  start 
time  of  each  cycle,  elevation  angle,  the  seven  reduced  values  of  CnS  (at 
nominal  altitudes  above  the  observatory  of  2.25,  3.75,  5.25,  7.5,  9.  i 5, 
12.75,  and  > 14.  5 km),  the  twenty  minute  mean  log -amplitude  variance  and 
the  standard  deviation  to  mean  ratio.  As  an  indication  of  the  relative  amount 
of  non- sta tiona r ity , the  ratio  of  the  number  of  valid  profiles  to  the  total 
collected  (yield)  is  included  in  Table  B-l.  Graphs  of  the  nightly  average  of 
the  valid  profiles  are  also  included  in  the  Appendix.  While,  in  fact,  the 
instrument  calculates  a value  appropriate  to  an  altitude  range,  the  results 
are  plotted  as  a function  of  the  peak  of  the  various  weighing  functions  used 
in  the  theory.  (12)  Because  of  the  small  zenith  angles  involved,  the  data 
has  not  been  corrected  for  this  effect. 

Several  interesting  properties  can  be  seen  in  the  results.  For  many 
of  the  profiles  lower  values  of  Cn2,  are  indicated  at  intermediate  levels. 
However,  these  levels  also  tend  to  fluctuate  more  rapidly  in  time  and  have 
relatively  greater  magnitude  variations  than  the  high  and  low  levels.  Not 
only  do  the  results  vary  substantially  from  night  to  night  but  large  changes 
also  occur  during  a single  night.  Noticeable  changes  are  seen  in  some 
cases  between  two  profiles  adjacent  in  time.  Examples  of  this  are  shown 
in  Figures  7 and  8. 

In  Figure  9,  the  average  profiles  for  each  month  in  which  data  was 
collected  are  shown.  Again,  they  are  plotted  as  a function  of  the  nominal 
altitudes.  The  number  of  profiles  averaged  for  each  curve  varies  as  a 
result  of  the  varying  amount  of  data  collected  during  each  month. 

3.  SCINTILLATION  DATA 

Measurements  of  the  aperture  averaged  log-amplitude  variances 
were  obtained  on  thirty  nights  during  the  experimental  period.  Table  B-l 
summarizes  these  data  giving  the  date,  experimental  time  period,  number 
of  samples  and  the  mean  and  variance  of  the  twenty  minute  averages  obtained 
on  each  occasion.  Because  most  of  the  data  was  collected  at  zenith  angles  of 
less  than  30°,  corrections  were  not  made  for  this  effect.  The  mean  value 
for  the  entire  21  8 member  ensemble  is  5.85  x 10*4  with  a range  of  (1.6" 
28.3)  x 10"4. 

'j'o  provide  a more  detailed  look  at  the  statistics  of  scintillation, 
the  complete  set  of  24  second  (as  opposed  to  the  20  minute  data  of  Table 
B-l)  log-amplitude  variances  taken  during  the  single  four  hour  run  of 
18  November  1975  has  been  analyzed  as  a single  ensemble.  The  resulting 
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Star  Sensor  Profiles  of  18  November  1 975.  Altitude  is  in 
meters  above  the  site  and  the  numbers  above  each  profile 
indicate  the  measured  value  of  rQ  (in  cm)  obtained  with  the 
beeing  Monitor  within  each  cycle. 
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Figure  9. 


Star  Sensor  Profiles  - Monthly  Average  for  Nov-Dec  (1),  April  (2), 
May  (3),  June  (4)  and  July  (5).  Altitude  is  in  meters  above  the  site 
The  dashed  line  is  from  the  Hufnagel  wind  correlated  model  with 
an  average  wind  speed  of  15  m/s. 


histogram  and  cumulative  probability  density  are  shown  in  Figures  10  and  11. 
This  ensemble  contains  534  points  and  has  a mean  and  root  variance  of  3.  76 
x 10-4  and  0.93  x 10“  , respectively.  The  smooth  curves  represent  a 
Gaussian  distribution  with  mean  and  variance  equal  to  the  empirical  values. 

4.  STAR  SENSOR  - SMALL  APERTURE  PHOTOMETER  COMPARATIVE  DATA 

In  order  to  obtain  direct  and  simultaneous  observations  of  the  log- 
amplitude  variance,  the  photometer  system  on  the  AMOS  b-29  telescope 
was  equipped  with  an  equivalent  36cm  aperture  to  correspond  to  that  of  the 
Star  Sensor.  A central  obscuration  was  not  used,  but  this  should  be  only 
a small  effect.  The  b-29  data  was  reduced  using  the  PDP-8  data  processing 
system.  Operations  were  carried  out  during  nine  nights.  On  two  occasions 
(26  May  and  9 June  1976)  synchronous  twenty-four  second  average  variances 
were  obtained.  The  rest  of  the  data  corresponds  to  10  min  averaged  b-29 
data  taken  within  the  basic  20  min  cycle  of  the  Star  Sensor. 

The  complete  data  set  is  given  in  Figure  12.  If  the  agreement  were 
perfect,  all  data  would  lie  on  the  solid  line.  The  dashed  lines  indicate  a 
factor  of  two  difference  in  the  measurements.  While  there  is  some  disper- 
sion in  the  data,  the  agreement  is  reasonably  good. 

One  effect  which  is  not  accounted  for  in  the  data  of  Figure  12  is  the 
contribution  due  to  noise.  To  investigate  this  effect,  b-29  data  was  taken 
at  a variety  of  bandwidths.  The  results  of  one  such  data  run  (taken  on 
26  May)  is  shown  in  Figure  13.  Assuming  the  noise  is  due  to  photoelectron 
fluctuations  and  that  the  bandwidth  of  significant  atmospheric  fluctuations 
is  below  1000  Hz  (and  stationary),  the  higher  bandwidth  data  should  be  des- 
cribable  by  a constant  plus  linear  frequency  term.  The  smooth  curve  in 
the  figure  represents  a least  squares  fit  of  this  model  to  the  three  wide 
bandwidth  data  points.  The  fit  is  very  good  (better  than  5%).  This  result 
implies  that  the  noise  contribution  to  the  b-29  data  is  approximately  2 x 10"4. 
Similar  results  were  also  obtained  on  the  night  of  9 June,  however,  there 
was  more  dispersion  in  the  data.  Subtracting  this  level  of  fluctuation  from 
the  b-29  data  would  place  most  of  the  points  to  the  left  of  the  solid  line  in 
Figure  12.  There  is,  of  course,  also  a noise  contribution  to  the  Star  Sensor 
data.  Assuming  this  to  be  of  order  10'4,  one-half  of  the  derived  b-29  noise 
was  subtracted  from  the  data  and  the  24  sec  average  data  replotted.  The 
results  are  given  in  Figure  14.  It  should  be  noted  that  the  "adjustment"  was 
derived  from  one  data  set,  but  fits  the  other  data  set  equally  well.  As  can 
be  seen  from  this  figure,  most  of  the  data  lies  quite  close  to  the  equal  magni- 
tude solid  line.  The  dashed  lines  represnet  a-i-10%  difference.  The  maximum 
difference  seen  is  ( + 27.8,  -27.  4)%. 

The  conclusion  that  can  be  drawn  from  this  analysis  is  that  the  Star 
Sensor  log -amplitude  variance  has  a noise  contribution  of  order  10*4.  An 
independent  set  of  noise  evaluation  tests  carried  out  on  the  Star  Sensor 
(and  reported  in  the  next  section)  indicate  that  for  the  operating  conditions 
of  these  experiments,  the  noise  level  is  no  larger  than  0.  5 x 10"4.  While 
there  is  a discrepancy  of  approximately  a factor  of  two,  the  combined 
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Star  Sensor  Log  - Amplitude  Variance  of  18  November  1°7 
Histogram.  The  smooth  line  is  for  a Gaussian  distribution 
with  the  empirical  mean  and  variance. 
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Figure  1 1.  Star  Sensor  Log  - Amplitude  Variance  of  18  November  1 975  - 
Cumulative  Probability.  The  straight  line  is  for  a Gaussian 
distribution  with  the  empirical  mean  and  variance 


/<I> 


STAR  SEN< 


66192  ADJUSTED  b-29  PHOTOMETER 


Figure  14.  Comparative  Log  - Amplitude  Variances  - Adjusted  Data 

Data  includes  24  sec  averages  taken  on  26  May  (o)  and  9 June  (•) 
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results  certainly  established  considerable  confidence  in  the  basic  operation 
of  the  Star  Sensor,  particularly  when  the  measured  log -amplitude  variance 
is  significantly  above  a value  of  (0.  5-1 . 0)  x 10“4. 

5.  NOISE  EVALUATION  TESTS 


To  provide  for  a direct  estimate  of  the  Star  Sensor  noise  characteristics 
a second  and  different  set  of  measurements  were  carried  out.  For  these 
experiments,  a low  voltage  light  bulb  was  connected  to  a regulated  DC  power 
supply  and  located  such  that  it  illuminated  a portion  of  the  diffuse,  white 
interior  surface  of  the  dome.  The  Star  Sensor  telescope  was  pointed  at  the 
illuminated  dome,  the  tracking  motor  was  turned  off  and  the  mount  was  locked 
in  position.  The  instrument  was  operated  with  photomultiplier  voltages  and 
outputs  typical  of  the  range  of  conditions  seen  during  stellar  operations. 

A laboratory  photometer  was  located  near  the  light  source  and  its  analog 
output  was  connected  to  the  PDP-8  data  processing  system  to  provide  a 
measure  of  the  illumination  fluctuations.  All  other  lights  were  turned  off 
and  personnel  were  excluded  from  the  dome  during  periods  of  data  collection. 
Data  was  collected  on  six  nights  over  a period  of  two  weeks  with  the  spatial- 
frequency  scanning  motor  both  on  and  off. 

The  results  of  these  measurements  are  shown  in  Figures  15  and  16. 

E(l)  is  the  value  of  the  total  aperature  log-amplitude  variance  while  E(2)  is 
the  log-amplitude  variance  associated  with  the  spatially  filtered  signal  (i.e.  , 
the  data  required  for  profile  reduction).  E(3)  is  the  average  photomultiplier 
output  (in  arbitrary  units)  and  the  voltage  given  is  that  applied  to  the  photo- 
multipliers. The  straight  lines  in  Figure  15  represent  a lea st- squa res  fit 
to  form 

E(2)  = A ebV  [ E ( 3 ) ] " 1 + C 

The  inverse  proportionality  to  E(3)  is  consistent  with  an  instrumental 
noise  dominated  by  photoelectron  fluctuations.  The  constant  term  is  assumed 
to  be  a contribution  due  to  the  light  source.  In  the  least-square  fit  several 
points  which  obviously  did  not  fall  close  to  the  straight  line  were  excluded. 

The  data  actually  plotted  has  been  reduced  by  the  empirically  derived  constant 
value.  As  can  be  seen,  the  fit  is  quite  good  except  for  a few  points.  The 
empirical  fitgave  A = 6.75  x 10*6  and  B = 0.  01037. 


The  assumed  source  contribution  (C)  was  found  to  be  in  the  range 
(1-3)  x 10“^  except  for  the  data  taken  on  25  August  which  was  of  order 
5.5  x 10"^.  The  lower  values  are  consistent  with  the  estimates  obtained 
from  the  photometer.  There  was  indications  that  due  to  unknown  causes,  the 
line  voltage  in  the  dome  was  low  and  variable  during  the  nights  of  2 3 and 
25  August  which  could  have  effected  the  light  source  and/or  the  Star  Sensor 
and  resulted  in  a higher  level  of  noise.  This  level  of  source  noise  could  also 
explain  the  higher  dispersion  in  the  data  at  low  voltage  and  high  E(3). 


The  straight  lines  in  Figure  16  also  represent  a least-square  fit  to 
the  model  given  in  Eq.  (7).  This  fit  gives  A=  2.47x10"-  and  B = 0.  01075. 
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TOTAL  LOG- AMPLITUDE  VARIANCE  : E(l) 
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Figure  15.  Star  Sensor  Noise  - Total  Log  - Amplitude  Variance.  Included 
is  motor  off  data  of  25  August  (•>.  26  August  (x),  27  August  ( + ), 
31  August  (o)  and  1 September  (Q)  as  well  as  motor  on  data  of 
23  August  (V),  31  August  (A)  and  1 September  (□).  The  straight 
lines  are  least-square  fits  of  Eq.  (7)  to  the  data  for  the  PM 
voltage  indicated.  The  three  points  marked  (NOAA)  are  the 
values  used  in  the  noise  compensation  subroutine  at  the  voltage 
indicated. 
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FILTERED  LOG-  AMPLITUDE  VARIAMC 


Figure  16.  Star  Sensor  Noise  - Filtered  Log  - Amplit 
as  Fig.  15. 
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The  empirical  va lue s of  C were  much  smaller  (below  10"^)  and  hence  were 
set  to  zero.  This  includes  the  data  of  25  August  which  showed  relatively 
higher  levels  in  the  values  of  E(l).  As  can  be  seen,  there  are  also  fewer 
points  well  off  the  curves  than  in  Figure  15.  All  of  these  observations  are 
very  consistent  with  the  interpretation  of  the  constant  factor  as  a light 
source  fluctuation.  The  reason  for  this  is  that  E(2)  is  derived  from  the 
difference  in  outputs  of  the  two  photomultipliers,  each  collecting  approx- 
imately one-half  of  the  light.  Hence,  any  overall  temporal  (but  spatially 
stationary)  fluctuation  in  the  illumination  should  cancel. 

Another  issue  is  the  spatial  uniformity  of  the  illumination.  If  this 
effect  were  present,  E(2)  might  indeed  be  measuring  a real  spatial  variation 
which  should  have  varied  as  the  spatial  frequency  was  changed.  No  evidence 
of  this  was  seen  as  is  indicated  by  the  data  taken  with  the  filter  drive  motor 
on. 


Experiments  similar  to  these  were  carried  out  by  N'OAA  before  the 
instrument  was  shipped  to  AMOS  and  a noise  cancellation  routine  was  pro- 
grammed into  the  software.  Only  a voltage  dependence  was  included.  The 
numerical  values  are  indicated  in  the  figures.  The  constants  (b)  are  quite  close 
(0.01103  for  E(2)  and  0.0126  for  E(2))  but  the  overall  level  of  noise  is  some- 
what higher  in  the  present  experiment.  While  there  is  no  definite  explanation 
for  this,  it  could  be  attributed  to  a variety  of  causes  including  different 
operating  environment,  aging,  etc.  However,  the  illumination  level  depend- 
ence found  in  the  present  experiments  was  expected  and  probably  should  be 
incorporated  into  the  noise  cancellation  routine. 

Complete  twenty  minute  cycles  were  also  run,  resulting  in  turbulent 
profiles  derived  in  this  experimental  configuration.  While  differences  are 
seen  in  profiles  taken  in  sequence,  insufficient  data  exist  to  obtain  a real 
estimate  of  the  errors  associated  with  the  Cn  values.  Subjective  impressions 
tend  to  support  a il  0%  accuracy  but  a more  complete  theoretical  and/or 
experimental  assessment  is  required  to  establish  any  firm  conclusion. 

A rough  estimate  of  the  range  of  reliable  operation  of  the  Star  Sensor 
can  also  be  obtained  from  these  results.  The  operations  carried  out  over  the 
past  year  were  surveyed  to  obtain  an  estimate  of  operating  conditions  as  a 
function  of  stellar  magnitude.  The  results  for  applied  PM  voltage  vs. 
magnitude  are  given  in  Figure  17.  The  operating  range  for  E ( 3 ) was  typically 
(600-900).  These  results  imply  that  a magnitude  3 star  would  require  an 
operating  voltage  of  900  V while  a magnitude  4 star  would  require  1000  V 
(maximum  allowed).  The  results  of  Figure  15  indicate  that  for  the  E(3) 
operating  range  the  noise  contribution  to  the  total  log -amplitude  variance  would 
be  (1.4-2. 1)  x 1 0"4  and  (2. 4-3. 6)  x 10 "4  for  a magnitude  3 and  4 star, 
respectively.  These  values  are  within  the  range  of  measured  atmospheric 
log-amplitude  variances.  While  the  data  reduction  software  does  subtract 
out  the  estimated  noise  contribution,  an  experimental  result  which  has  a 
spurious  contribution  of  magnitude  equal  to  that  of  the  signal  should  be 
viewed  with  caution.  This  is  particularly  true  of  the  Star  Sensor  because 
a complete  evaluation  of  the  effect  of  noise  on  the  profile  is  not  yet  available. 
Therefore,  for  operational  purposes,  the  threshold  for  useful  data  has  been 
set  at  stellar  magnitude  3 to  4. 
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V.  METEOROLOGICAL  DATA 


The  AMOS  Observatory  is  equipped  with  two  18-meter  towers  instru- 
mented with  routine  meteorological  instruments  and  microthermal  probes. 
Some  of  the  meteorological  data  is  necessary  for  the  reduction  of  micro- 
thermal  and  acoustic  sounder  data.  In  addition,  routine  measurements  are 
desirable  for  establishing  the  gross  meteorology  of  the  site.  While  some  of 
this  latter  data  hs  been  collected  and  will  be  in  the  future,  results  are  not 
reported  here . 

Mounted  on  each  tower  is  a triad  of  microthermal  probes  with  a 
nominal  one  meter  separation.  The  analog  outputs  from  these  probes  are 
connected  to  the  PDP-8  data  processing  system  which  calculates  a mean, 
variance  and  covariance  matrix  for  each  array.  The  latter  two  quantities 
are  required  for  the  calculation  of  and  Cn^  values.  For  all  data  reported 

the  averaging  time  was  set  at  approximately  10  min.  While  six  values  are 
calculated,  one  is  preferred.  The  selection  was  based  on  the  probe  pair 
whose  axis  was  closest  to  normal  to  the  wind  direction  on  the  windward 
tower.  The  first  requirement  is  based  on  geometrical  considerations  while 
the  second  is  to  avoid  a potentially  turbulent  wake  due  to  the  observatory 
domes.  However,  an  often  overriding  consideration  to  pair  selection  was 
an  obviously  broken  sensor.  If  this  occurred,  the  single  active  pair  remaining 
on  the  windward  tower  was  used. 

While  data  was  collected  on  a number  of  occasions,  only  that  obtained 
during  November  and  December  of  1975  is  reported.  There  were  several 
reasons  for  this.  One  is  that  this  period  represents  a major  data  collection 
effort  while  the  1976  measurements  were  more  routine.  More  importantly, 
was  the  failure  of  the  acoustic  sounder  to  provide  reduced  data.  Without 
the  additional  profiling  device,  a single  value  of  Cn^  at  ground  level  was 
judged  to  be  of  marginal  significance,  particularly  when  viewed  relative  to 
maintenance  and  manpower  requirements. 


Appendix  C summarizes  the  data.  Table  C-l  gives  the  dates,  time 
period,  number  of  samples  and  mean  value  and  range  of  the  reduced  Cn^ 
values.  The  graphs  show  the  temporal  behavior  of  the  reduced  data.  Foi^ 
the  total  data  set  of  143  points,  the  mean  value  of  Cn^  is  5.26  x 10"^  M"1-' 

As  can  be  seen  from  the  various  figures,  temporal  behavior  varied  from  a 
relatively  constant  value  over  a three  hour  period  (17  November)  to  fluctuations 
over  a similar  time  period  of  almost  a factor  of  10  (19  and  21  November). 

One  run  (14  November)  show  a tendency  towards  increasing  values  of  Cn*- 
while  one  run  (6  December)  shows  a tendency  towards  decreasing  values 


of  C. 


However,  the  existing  data  set  is  not  sufficient  to  establish  any 


firm  conclusions  regarding  these  types  of  effects.  When  compared  to 
direct  optical  data,  no  obvious  correlation  is  apparent. 
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Also  included  in  the  Appendix  is  gross  meteorological  information 
obtained  from  radiosonde  data  of  the  U.S.  Weather  Bureau  at  Hilo,  Hawaii 
a distance  of  approximately  100  mile 6 southeast  of  AMOS.  The  launch 
time  was  0200  HST  and  the  dates  indicate  the  mission  date  as  given  in 


Table  B-2.  The  existence  of  an  upper  atmospheric  disturbance  from 
17  November  to  23  November  is  indicated  by  the  presence  of  a jet  stream 
at  approximately  12,000  M (this  jet  is  usually  found  farther  north)  and  the 
replacement  of  the  usual  northeast  trade  winds  by  roughly  southwest  winds 
at  altitudes  below  12,000  M.  While  detailed  studies  of  the  correlation  of 
this  type  of  data  with  optical  data  has  not  been  carried  out,  some  comparisons 
have  been  made  and  are  reported  in  Section  VI. 
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VI.  DATA  COMPARISONS 


1 . GENERAL 

In  this  section  a number  of  comparative  studies  of  the  data  are  re- 
ported. Several  types  are  included.  The  first  is  a direct  comparison  of 
optical  parameters  with  theoretical  values  calculated  on  the  basis  of  the 
empirical  turbulence  profile  information.  Others  include  Seeing  Monitor- 
Star  Sensor  non- stationarity  and  variability,  experimental  vs  theory  and 
comparisons  with  other  experimental  data  reported  in  the  literature.  In 
general,  the  data  appears  to  be  consistent  in  some  aspects  but  with  poor 
agreement  between  the  profile  and  direct  optical  measurements. 

2.  TURBULENCE  PROFILE 

Turbulence  profiles  averaged  over  the  complete  data  runs  taken  on 
17,  18  and  21  November  1975  are  shown  in  Figures  18,  19  and  20.  The 
lowest  altitude  data  points  were  obtained  from  the  ground  based  microthermal 
sensors  (Appendix  C).  Data  from  37  m to  2.5  km  was  provided  by  NOAA 
from  a reduction  of  airborne  microthermal  data  they  collected  using  an 
instrumented  light  aircraft.  (H)  The  line  segments  from  1 km  to  24  km 
were  derived  from  Star  Sensor  profiles  (Appendix  B)  and  represent  the 
approximate  widths  of  the  weighing  functions  used  in  the  data  reduction.  ( 1 2) 
The  horizontal  scale  is  in  height  (meters)  above  the  observatory.  The  trop- 
opause  height  was  estimated  from  the  temperature  profiles  obtained  from  the 
U.  S.  W eather  Bureau  (Appendix  C). 

Values  of  rQ  and  o \Z  (point  aoerture)  were  calculated  from  these 
profiles  using  the  expressions  H),  (13), 
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13.  R.  S.  Lawrence  and  J.  W.  Strohbehn,  Proc.  of  IEEE  58,  1523(1970). 
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Figure  18  Empirical  Turbulence  Profile  of  17  November  1 975.  Altitude 
is  in  height  above  the  observatory.  Based  on  ground  level 
microthermal  data  (□),  airborne  microthermal  data  (+ ) and 
Star  Sensor  data  (i — — i).  Estimated  tropopause  height  of 
11  to  15.5  km. 
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The  calculation  was  carried  out  (at  RADC)  on  a Hewlett  Packard 
9810  desk  top  calculator  using  a 45  point  Riemann  sum  in  the  region  be- 
tween 15  m and  25  k.  VYhere  experimental  data  overlaps,  the  Star  Sensor 
results  were  used.  All  results  are  scaled  to  a wavelength  of  5000  A and 
are  given  in  Table  2. 

. 7 

Comparative  values  of  rQ  and  aj  (D  = 36  cm)  are  given  in 
Table  3.  Profile  values  of  ctj  2 (D  = 36  cm)  were  obtained  from  the  re- 
sults of  Table  2 by  multiplying  by  the  theoretical  aperture  averaging 
factor  of  0.  019.  04)  The  Seeing  Monitor  values  are  slightly  higher  than 
those  reported  in  Appendix  A.  The  reason  for  this  is  two-fold.  The  first 
is  that  all  values  were  corrected  for  zenith  angle  and  the  second  is  that 
each  value  of  rD  was  reduced  to  the  equivalent  value  of  the  Cn^  integral 
before  averaging  to  obtain  a more  direct  comparison  with  the  profile  in- 
formation. Table  3 also  includes  the  factor  by  which  the  rQ  profile  inte- 
gral would  have  to  be  multiplied  to  obtain  agreement. 

Comparisons  between  Star  Sensor  profiles  and  the  direct  measure- 
ment of  ai^  have  been  carried  out  by  NOAA.  ( The  results  are  given  in 
Figure  21.  These  two  values  are  independent  because  except  for  a minor 
normalization  proceedure,  the  measured  log -amplitude  variance  is  not  used 
in  the  profile  data  reduction. 

As  can  be  seen  from  these  results,  agreement  between  the  profile 
and  Seeing  Monitor  data  is  not  good.  A possible  explanation  for  this  is 
that  significant  amounts  of  turbulence  was  missed  by  the  profiling  instru- 
mentation. Assuming  this  to  be  the  case,  the  good  agreement  of  the  scin- 
tillation data  would  indicate  that  the  missing  turbulence  is  at  low  altitudes. 

It  should  be  noted  that  the  low  altitude  airborne  data  was  not  taken  directly 
over  the  observatory  but  over  a region  on  a lower  portion  of  the  mountain. 
Therefore,  the  data  obtained  may  not  be  representative  of  conditions  at  the 
observatory  because  the  relative  distance  to  the  ground  was  greater. 

The  required  increase  of  the  Cn^  integral  to  bring  the  two  values 
into  agreement  could,  for  example,  be  accounted  for  by  a turbulent  layer 
of  order  20-80  meter  thick  of  average  strength  Cn2  = 10*^  m'2/3.  These 
values  are  not  unreasonable.  If  such  a layer  existed  in  the  first  300  meters 
above  the  observatory,  it  could  have  been  seen  by  the  acoustic  sounder.  Un- 
fortunately, quantitative  reduction  of  this  data  was  not  successful. 

The  possibility  that  spatial  variability  of  turbulence  is  a potential 
cause  of  the  discrepancy  is  clearly  demonstrated  by  the  data  of  Figures 
22  to  24.  These  profiles  were  obtained  by  NOAA  from  airborne  data  they 


14.  H.  Yura,  Aerospace  Technical  Report  #TR0077(2756)- 1 , 1976. 

15.  R.S.  Lawrence,  private  communication. 
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TABLE  2.  OPTICAL  PARAMETERS  OBTAINED  FROM 
EMPIRICAL  PROFILES 


Date 

Correlation  Scale 
rQ  (cm) 

Log -Amplitude  Variance 
(point  aperture) 

17  November 

19.  9 

0.  028 

1 8 November 

19.  2 

0.  012 

21  November 

12.  8 

0.  029 
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Figure  21.  Comparison  of  Measured  and  Profile  Derived  Lop  - Amplitude 
Variances.  Only  valid  profiles  are  used.  Big  dots  are  nightly 
averages  and  the  small  dots  are  for  single  2h  min  cycles 


West  Sea  Level,  (km)  East 


Figure  22.  Two  Dimensional  Turbulence  Contours  of  17  November  1075 
Based  on  daytime  microthermal  data.  The  observatory  is 
approximately  0.  5 km  west  of  the  peak.  Winds  are  from  the 
trade  direction. 


Figure  24.  Two  Dimensional  Turbulence  Contours  of  20  November  1975 
Same  as  Fig . 22 . 


collected  during  the  daytime.  If  nighttime  conditions  are  similar,  then 
it  is  quite  possible  that  the  aircraft  sensor  was  sampling  substantially 
different  turbulence  than  the  optical  measurements.  In  these  profiles,  the 
observatory  is  located  slightly  to  the  west  of  the  peak  (i  0.  5 km)  and  the 
winds  are  from  the  predominant  trade  direction  (northeast). 

3.  SEEING  MONITOR  - STAR  SENSOR 


Several  comparisons  can  also  be  made  between  Seeing  Monitor  and 
Star  Sensor  data.  One  such  comparison  is  given  in  Figures  7 and  8.  The 
numbers  given  at  the  top  of  each  profile  are  the  optical  values  of  r0  ob- 
tained during  a 10  min  data  collection  cycle  within  the  20  min  Star  Sensor 
cycle.  There  appears  to  be  a rough  correlation  between  these  two  results 
(i.  e.  , rQ  larger  for  profiles  with  smaller  area),  at  least  for  the  run  of 
18  November.  Of  course,  a strong  correlation  would  be  expected  only  if 
the  high  altitude  portion  of  the  profile  dominates  and/or  the  turbulence  levels 
under  1 km  are  relatively  constant.  The  former  assumption  is  generally  not 
valid  and,  in  fact,  the  data  supports  this  conclusion.  The  average  profiles 
given  in  Figure  9 yields  a value  for  the  integral  of  Cn  for  altitudes  above 
1 km  of  the  order  of  2 x 10-13  m1/3.  Using  the  theoretical  expression  for  ru 
and  the  average  value  of  9.  6 cm  at  5000  K yields  a value  for  the  integral  of 
Cn2  over  the  entire  atmospheric  path  of  approximately  7 x 10'13  m1/  3. 
Therefore,  based  on  these  averages,  turbulence  above  1 km  accounts  for 
roughly  35%  of  the  total.  Furthermore,  based  on  the  airborne  microthermal 
data,  the  region  between  300  m and  1 km  contributes  roughly  10"1:>  m4  1 ' - 
or  20%  of  the  total.  These  values  lead  to  the  conclusion  that  approximately 
4 5%  of  the  total  contribution  to  the  integrated  turbulence  is  associated  with 
low  levels  (under  300  m). 


Correlation  also  exist  in  non-stationary  effects  seen  in  the  two  data 
sets.  For  example,  the  Star  Sensor  yields  for  6 and  8 December  were  82 
and  30%,  respectively  (see  Figure  4 for  Seeing  Monitor  data).  On  a cycle 
by  cycle  basis,  the  standard  deviation  to  mean  ratio  for  the  Seeing  Monitor 
output  voltage  had  an  average  of  0.  08  and  a range  of  (0.  05-0.  16)  on  6 De - 
c embe r.  For  8 December,  the  values  were  0.  17  and  (0.  08-0.  21).  The 
comparable  Star  Sensor  data  is  the  twenty  minute  standard -deviation  to 
mean  ratio  of  the  log -amplitude  variance.  The  average  and  range  of  tnis 
pa rametcr  wa s 0.  16  (0.  11-0.  21)  and  0.27  (0.  13-0.46)  for  6 and  8 Dec emb e r , 
respectively.  Hence  both  sets  of  data  display  higher  levels  of  variability 
on  8 December. 


Similar  trends  are  sometimes  apparent  in  the  data.  For  example, 
the  Star  Sensor  log -amplitude  variance,  except  for  one  point,  decreased  in 
a monotonic  fashion  during  the  run  of  12  July  indicating  improving  conditions 
in  agreement  with  the  Seeing  Monitor  data  of  Figure  5.  In  contrast,  no  ob- 
vious trend  was  seen  in  the  log -amplitude  variance  during  the  run  of  8 lulv. 
Because  scintillation  tends  to  be  dominated  by  higher  altitude  turbulent  < 
while  seeing  (rD)  is  usually  dominated  by  lower  altitudes,  the  two  sets  ol 
data  should  not  necessarily  show  the  same  trends. 
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4.  THEORETICAL  COMPARISONS 


i ) 4 i 

H.  Yura'  has  carried  out  a number  of  theoretical  calculations 
with  regard  to  these  measurements.  The  theory  uses  a combination  of 
the  Hufnagel  wind  correlated  model(lk)  for  high  altitudes  and  a low  alti- 
tude model  based  on  the  work  of  BuftonO?)  and  Koprov  and  Tsuang.  ( 1 8 I 
The  results  of  these  calculations  are  given  in  Table  4.  The  input  to  the 
theoretical  model  is  the  average  wind  velocity  which  was  derived  from 
the  meteorological  data  of  Appendix  C.  The  empirical  values  for  this 
pa  ra meter  are  21.  3 m/s,  20.  5 m/ s,  and  30.  8 m/s  for  17,  18  and  21, 
November,  respectively.  As  can  be  seen,  the  theoretical  to  results  are 
higher  than  the  experimental  data  in  all  cases.  It  is  interesting  to  note 
that  for  rQ,  the  theoretical  values  are  in  better  agreement  than  the  values 
derived  from  the  empirical  profile  (Table  3). 

The  Hufnagel  model  can  also  be  compared  to  the  Star  Sensor  high 
altitude  profile  data.  In  Figure  9,  the  dashed  curve  is  based  on  the  model 
assuming  an  altitude  averaged  velocity  of  15  m/s.  Several  observations 
are  obvious.  All  of  the  experimental  results  have  a reasonable  shape.  In 
addition,  all  show  a tendency  to  flatten  out  at  intermediate  altitudes  and 
then  fall  off  rapidly.  Only  two  of  the  five  averaged  results  show  an  increase 
in  turbulence  strength  in  the  vicinity  of  (6-9  km)  above  the  site  (9-12  km 
above  sea  level).  In  general,  the  data  appears  to  support  the  existence  of 
a turbulent  tropopause  but  one  which  is  not  always  strong.  The  data  also 
indicates  relatively  higher  levels  of  turbulence  at  lower  altitudes  than  the 
theoretical  mode.  This  could  possibly  be  due  to  the  application  of  the  model 
outside  its  range  of  stated  validity  (ground  +3000  M to  24,  000  M ).  A rather 
dramatic  result  is  the  close  agreement  between  the  model  and  the  averaged 
data  for  July  1976.  Because  a detailed  study  of  the  available  meteorological 
data  on  upper  atmospheric  winds  has  not  yet  been  carried  out,  great  signif- 
icance cannot  be  attached  to  this  agreement. 

5.  EXPERIMENTAL  COMPARISONS 

Scintillation  data  has  been  collected  for  a considerable  period  of 
time  at  a number  of  different  sites.  A number  of  average  properties 
associated  with  this  data  base  are  given  in  Ref.  16.  Several  can  be  con  - 
pared  with  the  data  reported  here.  Results  are  given  in  terms  of  the 
parameter  S,  defined  as  ap^/  < p)>  ^ where  ap^  is  the  variance  and  <p> 
is  the  mean  value  of  the  in+egrated  irradiance  seen  by  a 10  cm  aperture 
viewing  near  zenith  stars.  In  particular,  over  a year  long  period  log  S 
is  approximately  a Gaussian  random  variable  with  a = log  2.  1.  The  year 
round  median  value  of  S is  about  0.  06  for  all  sites.  The  complete  set  of 

16.  R.  E.  Hufnagel,  OSA  Topical  Meeting  on  Optical  Propagation  Through 

Turbulence  (Boulder,  Colorado,  1974),  paper  WA1. 

17.  J.  L.  Bufton,  Appl.  Opt.  J_2,  1785  ( 1973). 

18.  V.  M.  Koprov  and  L.  R.  Tsuang,  Atmos,  and  Oceanic.  Phys. 

(1966). 
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TABLE  4.  THEORETICAL  COMPARLSONS 
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twenty  minute  log -amplitude  variances  given  in  Appendix  B has  a mean 
median  and  variance  of  approximately  5.  9 x 10-4,  5,  6 x 10*4  and  1 9 x'  10*7 
respective^.  However,  these  values  are  for  a 36  cm  aperture.  Assuming 
he  approximate  large  aperture  theoretical  scaling  (D’7/3)  and  the  usual 
uaussian  relationship  between  irradiance  and  log  -amplitude! 1 3)  yields  a 
median  value  for  S of  0.  04  and  a root  variance  for  log  S of  approximately 

Aunc  * number  of  other  measurements  of  rQ  have  been  carried  out.  At 
AMOS,  techniques  implemented  have  included  Speckle  Interferometry  ( 1 9 1 
a deferential  Hartmann  method,  (20)  and  a shearing  interferometer.  (21) 
ft?  a ttS  3 i°  co^ectec3  at  Ihe  Mauna  Kea  Observatory  (altitude  ^ 14  000 

of  AMOS  1Sn  ° fHawaii-  a distance  of  approximately  75  miles  southeast 
01  AMOS.  The  instrument  used  was  a coherent  interferometer.  ! 22 ) AU  of 
these  results  are  given  in  Table  5.  The  reported  shearing  interferometer 
ta  is  limited  to  a single  plot  of  the  long  exposure  MTF.  The  shape  varies 
somewhat  from  the  theoretical  expectation.  The  value  of  rQ  in  the^able 
presents  an  approximate  fit  of  the  theory  to  the  data.  The  speckle  data 

similar  to*  hose°n'StaV°?trity  effCCtS  and  aPProximate  Gaussian  statistics 
similar  to  those  reported  here. 


19. 


20. 

21. 

22. 


A.  Schneide rman  and  D.  Karo,  OSA  Topical  Meeting  on  Speckle 
Phenomena  (Pacific  Grove,  California,  1976),  paper  ThC6  and 
private  communication. 

M.  G.  Miller  and  P.  F.  Kellen,  AAS/SAO/OSA/SPSE  Topical  Meeting 
on  Imaging  in  Astronomy  (Cambridge,  Mass.  , 1975),  paper  WB  3 
D.  Kelsall,  J.  Opt.  Soc.  Am.  6_3,  1472  (1973). 

fl975  )DaintY  and  R* J’  Scaddan’  Mon-  Not.  R.  Astr.  Soc.  170,  519 
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TABLE  5.  EXPERIMENTAL  COMPARISONS  OF  tq 


T echnique 

Seeing  Monitor 

Speckle  Interferometry 
(Schneiderman  and  Karo) 

June  1975 

J un  e 1976 


Mean  (cm)  Range  (cm) 

9.6  5.3-17.8 

4.  1-6. 4 
10 


Differential  Hartmann, 

August  1974 
(Miller  and  Kellen) 

Coherent  Interferometer, 

June  1974 

(Dainty  and  Scadden -Mauna  Kea) 

Shearing  Interferometer,  1972 
(Kelsall ) 


3.  0-18.  0 


4.  0-18.  5 


3 
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VII.  R TAM  AMOS  TESTS 


1.  GENERAL 

The  Real  Time  Atmospheric  Measurement  System  (KTAM)  underwent 
AMOS  testing  during  the  period  16  August  to  3 September,  1976.  Initially,  it 
was  set  up  on  an  optical  bench  and  its  performance  checked  vising  a Ile-Ne 
laser  source.  Results  were  similar  to  those  seen  in  the  laboratory.  After 
installation  on  the  1. 6 M telescope,  a series  of  diagnostic  tests  and  adjust- 
ments were  carried  out.  Successful  operation  of  the  instrument  was  es- 
tablished. However,  several  problems  were  found  but  not  eliminated.  The 
first  of  these  was  a very  narrow  field  of  view  (less  than  fovir  arc  seconds) 
which,  if  not  widened,  will  be  a severe  restriction  to  any  future  operations, 
requiring  very  good  tracking  and  low  wind  velocity.  The  second  was  a low 
signal  to  noise  ratio  which  has  significant  impact  on  the  ability  of  the  R TAM 
to  produce  a detailed  characterization  of  the  OTF. 

2.  MTF  PERFORMANCE 

In  general,  the  MTF's  seen  were  consistent  with  expectation,  i.  e.  , 
a central  peak  which  fell  off  at  low  spatial  frequency  to  a low  signal  level 
which  included  large  variations.  W These  variations  are  due  both  to  (lie 
atmosphere  and  system  noise.  Examples  are  shown  in  Figvire  25.  The 
width  of  the  peak  can  be  used  to  obtain  a rovigh  estimate  of  seeing  conditions 
(i.  e.  , rQ),  while  the  relative  levels  of  the  zero  and  high  spatial  frequency 
portion  of  the  MTF  can  be  used  to  obtain  an  estimate  of  the  signal  to  noise 
ratio. 

Seeing  conditions  wore  determined  by  measuring  the  width:  of  the 
MTF  at  50%  of  its  initial  value.  A fit  of  this  frequency  to  theory  yields 
the  desired  estimate  of  r0.  Because  the  modulus  of  the  OTF  is  measured, 
the  appropriate  model  is  the  sho  rt  -exposure  averaged  MTF.  (4)  The  data 
analyzed  in  this  fashion  are  short  exposures  (single  realizations)  which  ar< 
corrupted  with  noise  and  hence  these  results  are  only  approximate.  The 
estimated  values  of  rQ  were  of  order  15  cm  on  3 1 August  and  in  the  range 
h to  12  i m on  2 September.  These  results  are  reasonable,  based  on 
previous  experience. 

l or  the  purpose  of  carrying  out  a detailed  characterization  of  the 
MTF,  an  important  consideration  is  the  inherent  signal  to  noise  ratio  of 
the  MTF'  at  high  spatial  frequencies  (i.  e.  , greater  than  the  atmospheric 
cutoff  fy\  ~ rQ).  If  this  quantity  is  too  small,  useful  information  will  be 
obtained  at  low  spatial  frequencies  only.  It  is  important  to  note  that  this 
high-frequency  (S/N)  is  different  than  the  (S/N)  considered  in  various 
HTAM  documents  which  is  more  characteristic  of  zero-frequency  behavior 
Provided  various  statistical  assumptions  are  made,  both  can  lie  related  to 
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VEGA  - 31  AUG 
Y - MTF 

MIN.  BANDWIDTH 


VEGA  - 25  SEP! 

Y - MTF 

MAX.  BANDWIDTH 


m in  imum 


the  ratio  of  the  average  MTF  to  the  average  noise  level  at  zero  spatial  fre- 
quency. In  the  current  case,  this  latter  quantity  can  be  estimated  from  the 
data  because  the  high-frequency  behavior  is  dominated  by  noise.  Theory^®) 
predicts  a high-frequency  (S/N)  for  a single  realization  of  order  one  due  to 
atmospheric  fluctuations  only.  A theoretical  estimate(23)  (assuming  an 
rc  of  1 0 cm)  of  the  multiplicative  factor  by  which  this  ideal  value  is  reduced 
by  noise  as  a function  of  the  zero-frequency  MTF  to  average -noise  ratio 
is  given  in  Figure  26.  Also  included  in  the  figure  is  the  estimated  range  of 
operation  seen  on  31  August  and  2 September  as  well  as  the  theoretically 
predicted  range  of  operation  against  a first  magnitude  star.  The  former 
was  established  by  analyzing  data  taken  against  several  bright  stars  while 
the  latter  is  based  on  theoretical  estimates  of  system  characteristics.  (24) 

Two  other  quantities  are  of  interest  in  determining  the  data  processing 
requirements.  The  first  of  these  is  the  number  of  frames  required  to  obtain 
a given  accuracy  in  the  final  average.  Assuming  that  independent  realiza- 
tions are  averaged  the  (S/N)  is  proportional  to  the  square  root  of  the  number 
of  frames.  The  second  is  the  fraction  of  the  total  average  high-frequency 
signal  which  is  due  to  the  MTF.  If  this  fraction  is  small,  processinc  be- 
comes more  difficult  due  to  the  existence  of  a large  spurious  signal.  These 
three  characteristics  of  the  data  are  tabulated  m Table  6 for  both  the  actual 
performance  seen  and  the  theoretically  estimated  performance. 

These  results  indicate  that  performing  at  its  current  level,  the  RTAM 
would  provide  data  which  would  be  extremely  difficult  to  process  for  useful 
MTF  information  at  high-spatial  frequencies.  This  is  particularly  true  at 
the  lower  end  of  the  operational  range  where  several  thousand  frames  of 
data  would  be  required  and  still  result  in  a high  frequency  sicnal  which  would 
be  more  than  90rc  spurious.  If  the  theoretically  predicted  performance  were 
to  be  achieved,  processing  should  yield  useful  results  for  brighter  stars 
under  reasonable  conditions. 

3.  PHASE  TRANSFER  FUNCTION 

The  situation  regarding  the  phase  outputs  is  more  difficult  to  analyze, 
tthile  the  RTAM  did  generate  signals  in  these  channels,  their  interpretation 
is  not  obvious  for  a number  of  reasons.  Examples  of  the  PTF  taken  against 
Vega  on  2 September  are  shown  in  Figure  27.  Maximum  excursions  seen 
were  of  order  + 1 0 it  radians  with  rms  values  of  order  2 r radians  These 
values  do  not  conflict  with  any  fundamental  understanding.  However,  a 
number  of  non-atmospheric  effects  are  probably  contributing  including  a 
systematic  error,  tracking  induced  image  wander,  noise  and  loss  of 
phase-lock  in  the  detection  process.  The  systematic  error  was  observed 
in  the  laboratory.  While  its  presence  complicates  processing,  it  can  be 
eliminated  (if  it  is  constant  in  time)  by  subtracting  a known  phase  las  a 


l.  3 . Based  on  an  analysis  si  m i 1 a r to  th  at  given  by  F.  Roddi  er  (reference 

20,  paper  ThC6),  A.  Schneiderman  and  D.  Karo  (reference  1 1,  pg  70' 
and  M.  G.  Miller  (reference  19,  paper  ThC4). 

24.  A.  MacGovern,  private  communication 
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REDUCTION  IN  H.  F.  (S/N) 


[mtf/noise]  at  zero-frequency 


Figure  26.  R TAM -Reduction  in  High-Frequency  (S/N) 


VEGA  - 2 SEPT. 

Y - PHASE 
MAX.  BANDWIDTH 


VEGA  - 2 SEPT. 
Y - PHASE 
MIN  BANDWIDTH 


• ■ 


G6I96 


RTAM-PTF  Data.  Source  was  Vega  (m  ^ 0)  on  2 September 
t'ppor  trace  is  with  maximum  bandwidth  and  lower  trace  is 
with  minimum  bandwidth. 


function  of  frequency)  from  each  realization  Tracking  induced  im.ipr 
wander  would  be  somewhat  more  complicated  to  handle  because  it  is  dy- 
namic. However,  the  dynamics  of  the  mount  are  such  that  little  motion  i in 
occur  with  a bandwidth  greater  than  10  Hr  Therefore,  if  subsets  of  data 
corresponding  to  total  collection  times  of  less  than  0.  1 sec  were  normalized 
to  eliminate  the  least-square  fit  of  the  average  phase  to  a linear  phase 
taper  plus  the  (assumed)  known  systematic  error,  much  of  this  telescope 
wander  should  be  eliminated.  If  this  procedure  did  not  prove  to  be  adequate, 
it  would  also  be  possible  (in  principle)  to  derive  the  required  wander  effect 
from  a time  - synchronized  recording  of  the  mount  shaft  angle  encoder  error 
s ignals . 

The  effect  of  system  noise  is  to  add  a statistical  fluctuation  to  the 
signal  which  is  already  random  due  to  the  atmosphere.  Because  one 
interesting  question  regarding  phase  is  the  type  of  statistics  it  obeys, 
noise  will  act  as  a spurious  signal  which,  if  large  enough,  can  mask 
the  atmospheric  effect.  Increasing  the  number  of  frames  reduced  will  not 
minimize  this  effect,  but  only  refine  the  estimates  of  the  combined  statis- 
tics. Hence  if  useful  information  is  to  be  obtained,  either  the  noise  must 
be  small  or  a set  of  noise  evaluation  experiments  and/or  theory  (as  a 
function  of  operational  parameters)  must  be  carried  out  so  that  this  effect 
can  be  deconvolved  from  the  combined  data.  The  amount  of  noise  seen  in 
the  MTF  would  suggest  that  the  former  is  an  unlikely  situation. 

The  most  difficult  problem  to  deal  with  is  the  loss  of  phase -lock  in 
the  detection  process.  This  occurs  whenever  the  MTF  drops  below  some 
critical  level  which  is  unknown  at  this  time.  The  result  of  this  phase-lock 
loss  is  that  a random  and  unknown  amount  of  phase  shift  is  added  to  the 
output.  In  fact  if  the  MTF  does  not  rise  above  this  threshold,  the  phase 
output  will  apparently  be  due  entirely  to  noise.  Because  the  effect  of  the 
atmosphere  is  to  depress  the  MTF  at  all  frequencies  but  the  very  lowest, 
the  possibility  exists  that  the  only  real  atmospheric  phase  information 
which  can  be  obtained  from  the  RTAM  will  be  at  spatial  frequencies  below 
the  seeing  induced  cutoff.  If  this  is  indeed  the  use,  the  phase  outputs  will 
be  of  little  utility. 
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APPENDIX  A 
SEEING  MONITOR  DATA 
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Notes 


All  data  points  correspond  to  the  average  value  of  a 1 350  member 
ensemble  of  short  exposures  (one  millisec)  taken  over  a period  of 
approximately  8.  25  min. 

All  data  except  that  of  12  and  13  May  were  taken  with  a 1200  R 
filter  centered  at  6200  R.  The  data  of  12  and  13  May  was  not 
spectrally  filtered.  Theoretical  scaling  was  used  to  reduce  to 
the  effective  wavelength  of  5000  R. 

The  data  of  9-13  July  has  been  corrected  for  zenith  angle  depend- 
ence using  the  theoretical  scaling.  The  rest  of  the  data  was  not 
corrected  because  small  zenith  angles  (less  than  30°)  were  used 
and  hence  the  correction  is  small. 

The  times  given  are  approximate  (±5  min)  local  Hawaiian  Standard 
Time.  The  dates  given  are  those  of  the  starting  time. 


60 


TABLE  A - 1 . SEEING  MONITOR  DATA  - SUMMARY 


Date 

Time  Period 
HST 

No.  of 
Points 

ro  <§ 
Mean 

5000  K : cm 
Range 

11  Nov  197  5 

2240 

- 0030 

6 

9.  6 

7.9  - 10.7 

12 

1950 

- 2240 

9 

9.  9 

8.  5 - 1 1.  1 

14 

2110 

- 2400 

9 

11.7 

10.8  - 12.  9 

1 5 

1930 

- 2230 

10 

11.4 

10.  2 - 13.  0 

17 

1950 

- 2240 

9 

7.  7 

6.  8 - 9. 8 

18 

2050 

- 2350 

9 

11.3 

9.7  - 13.  2 

19 

2020 

- 2310 

9 

8.  7 

6.9  - 10.5 

21 

1930 

- 2220 

8 

8.  6 

6.6  - 12.4 

6 Dec  197  5 

0100 

- 0350 

9 

9.4 

9.  0 - 9.  7 

8 

0200 

- 0450 

9 

14.  1 

10.  0 - 17.  8 

12  May  1976 

2250 

- 0040 

6 

10.  7 

10.1  - 11.4 

I 3 

2120 

- 2300 

6 

12.  1 

10. 3 - 13.1 

27 

2230 

- 2320 

3 

5.  6 

5.  3 - 5.9 

1 0 June  1976 

2110 

- 2220 

4 

7.  0 

5.  3 - 8.  1 

1 8 

2120 

- 2400 

10 

7.8 

7.1  -8.7 

21 

2030 

- 2230 

12 

7.  6 

7.9  - 9.0 

24 

2150 

- 2350 

12 

10.  5 

9.1  - 11.3 

29 

2250 

- 0050 

12 

12.  2 

11.6  - 13.  3 

30 

2120 

- 2320 

12 

6.  0 

5.  5 - 6. 9 

6 July  1976 

2300 

- 0100 

12 

7.  6 

6.  2 - 8. 4 

8 

2330 

- 0130 

12 

9.  6 

8.1  - 11.6 

9 

2100 

- 0210 

9 

8.  2 

7.1  -9.0 

1 2 

0150 

- 04  50 

1 8 

11.1 

7.6  - 14.4 

1 3 

0050 

- 0420 

13 

8.  7 

8.  0 - 9. 6 
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Figure  A - 1 . Seeing  Monitor  Data  - 1 1 , 12,  1 4 and  1 5 November  1975  . 

The  horizontal  bars  indicate  the  averaging  period 
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Seeing  Monitor  Data  - 6,  8,  9,  12  and  13  July  1976.  Same 
as  A - 1 . 
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Notes 


The  basic  data  obtained  corresponds  to  ensembles  of  approximately 
40  twenty-four  second  averages  taken  during  a twenty  minute  time 
period.  These  data  are  then  reduced  for  the  turbulence  profile  and 
20  min.  averaged  a j’. 

The  reported  oj  2 is  the  mean  and  range  of  the  20  min.  averages 
taken  on  each  night. 

The  times  given  are  approximate  local  Hawaiian  Standard  Time. 

The  dates  given  are  those  appropriate  for  the  beginning  of  each 
night's  data  collection  run. 

The  data  has  not  been  corrected  for  zenith  angle.  Typically  this 
angle  was  less  than  30°  so  that  the  correction  should  be  small. 

No  spectral  filtering  was  used. 

Yield  refers  to  the  number  of  valid  profiles  divided  by  the  total 
number  or  profiles  collected  during  an  evening.  To  be  valid,  a 
profile  must  have  a standard  deviation  to  mean  ratio  of  the  twenty - 
four  second  log-amplitude  variances  of  0.  20  or  less.  Values 
larger  than  this  indicate  a high  level  of  non  - stat  iona  r ity  in  the 
atmos-phere  during  the  twenty  minute  collection  cycle.  Only  valid 
profiles  are  reported. 

No  valid  profiles  were  collected  on  11  and  12  November  due  to 
mechanical  and  electrical  failures  in  the  device.  Seven  of  the 
ten  profiles  taken  on  19  November  were  invalid  due  to  saturation 
of  the  photo-multiplier  tubes. 

The  seven  value  of  Cn  (i),  i = 1.  . . 7,  are  for  the  nominal  altitudes 
(above  the  site)  of  2.  25,  3.75,  5.25,  7.5,  9.75.  1 2.  75  and  > 14.  5 
km,  respectively.  SD/M  is  the  standard  deviation  to  mean  ratio  of 
the  ensemble  of  twenty-four  second  variances. 
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TABLE  STAR  SENSOR  DATA  - SUMMARY 


Date 

T ime 

Pe  riod 

Yield 

No.  of 

L. 

°1 

(x  10 

) 

HST 

Sample  s 

Range 

Mean 

11  Nov  1975 

2200 

- 

0040 

0/8 

- 

- 

- 

12 

1920 

- 

2220 

0/9 

2 

5.  50 

- 

5.  54 

5.  52 

14 

2010 

- 
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1 5 

2. 78 

- 

6.  22 

4.  16 

15 

1930 

- 

2310 

10/11 

13 

1.55 

- 

3.  98 

2.  68 

17 
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- 
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4/1  1 

12 
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18.  2 

7.92 
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1 1/12 

12 
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3.  76 

19 

1940 
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0/10 
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28.  3 

1 8.  2 

21 

2010 
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2250 

5/7 
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5.  17 
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8.  21 

6.00 

6 Dec  1975 

0020 

- 

0400 

9/1  1 

12 

4.  51 

- 

10.  4 

6.  88 

8 
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- 
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3/10 

10 

2.  1 3 

- 

5.  71 

3. 24 

1 6 April  1976 
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- 
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7/12 

12 
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- 

12.  8 

8.  1 3 

23 
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- 
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1 2/16 
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- 
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4.45 

29 
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4.83 

5 May  1976 
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6 

3.  54 

- 

8.  25 

6.  08 

12 
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- 

5.  22 
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2.  1 2 

26 
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- 
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1 
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6.  38 
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- 
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10 
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- 
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- 
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- 
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e, 
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21 . 3 

1 7.  8t 
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3.  9 1 

1 3 
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1 9 June  '76  2220  65  66.9  43. 
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Figure  B-l.  Star  Sensor  Nightly  Averaged  Profiles  - November  1975. 

Altitude  is  in  height  above  the  observatory.  Each  profile 
is  an  average  of  all  valid  data  collected  on  a specific  night. 
The  width  of  the  weighing  functions  are  not  shown. 
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Figure  B-4.  Star  Sensor  Nightly  Averaged  Profiles  - May  1976.  Same 
as  R - 1 . 
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Figure  B-6.  Star  Sensor  Nightly  Averaged  Profiles  - July  1976.  Same 
as  B - 1 . 
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APPENDIX  C 

METEOROLOGICAL  DATA 
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Note  s 


1.  All  data  points  correspond  to  the  average  value  of  a 1 350  member 
ensemble  taken  over  a period  of  approximately  8.25  min. 

2.  The  times  given  are  approximate  (+  5 min.  ) local  Hawaiian  Stan- 
dard Time.  The  dates  given  are  those  of  the  starting  time. 

3.  Only  one  of  the  six  calculated  values  of  Cn^  is  reported.  This 
value  was  selected  on  the  basis  of  wind  direction  and  proper 
operation  of  the  probes.  The  upwind  probe  pair  closest  to  nor- 
mal to  the  wind  direction  is  the  preferred  set. 

4.  Radiosonde  data  is  from  the  USWB  station  at  Hilo,  Hawaii,  ap- 
proximately 110  miles  south-east  of  AMOS.  Launch  time  was 
0200  HST.  Data  reported  corresponds  to  the  night  of  each 
data  run. 
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TABLE  C-l.  MICROTHERMAL  DATA  -SUMMAI 
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Figure  C - 3. 
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